Off-critical Casimir effect in Ising slabs for symmetric 
boundary conditions in spatial dimension d = 3 
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Extended de Gennes-Fisher (EdGF) local-functional method has been applied to the thermody- 
namic Casimir effect away from the critical point for systems in the Ising universality class confined 
between parallel plane plates with symmetric boundary conditions (denoted (ab) — (++)). Re- 
sults on the universal scaling functions of the Casimir force W++(y) (y is a temperature-dependent 
scaling variable) and Gibbs adsorption G(y) are presented in spatial dimension d = 3. Also, the 
mean-field form of the universal scaling function of the Gibbs adsorption G(y) is derived within the 
local functional theory. Asymptotic behavior of W++(y) for large values of the scaling variable y is 
analyzed in general dimension d. 

PACS numbers: 64.60.F-, 64.60.fd, 05.70. Jk, 68.15.+e, 68.35.Rh 



The Casimir effect in quantum or statistical physics 
refers to long-range forces that emerge due to a confine- 
ment on fluctuations. In statistical physics these fluctu- 
ations are in the order parameter of a thermodynamic 
system at or near the critical point, as predicted in 1978 
[l|. The Casimir force (CF) depends on the nature of 
the confined system as well as the boundary conditions 
(b.c.) and the geometrical form of the confinement 
Much theoretical work has examined various surface uni- 
versality classes for Ising systems and classical fluids ei- 
ther exactly at criticality or away from it [H . Symmetry 
breaking b.c. (defined below) are of particular interest for 
experiments with critical binary liquid mixtures. Appre- 
ciable agreement pertaining to these systems has been 
achieved between theory H,H, [1[ and recent experiments 
on the CF at criticality in complete wetting films of bi- 
nary /fluorocarbon mixture near liquid vapor coexistence 
@, with the mean value of the universal Casimir ampli- 
tude (a measure of CF at the bulk critical point, defined 
below) most closely corresponding to earlier prediction 
of the local- functional theory [4j, while at the same time 
encompassing other theoretical/simulation estimates 
Exact results on the full-temperature dependence of the 
CF are available in spatial dimensions d = 2 [7[ and d > 4 
(mean field theory) 0, Q for both cases of symmetry 
breaking b.c. 

Although significant theoretical effort has focused on 
the universal scaling functions of the off-critical Casimir 
effect, knowledge of them is still somewhat incomplete for 
spatial dimension d = 3, even for the relevant Ising uni- 
versality class. Pertinent results in this case refer to films 
with periodic b.c, studied via Monte Carlo (MC) simula- 
tions 0] or by the field-theoretic approach for Dirichlet, 
Neumann and periodic b.c. [Toj |. besides recent MC re- 
sults that now include symmetry breaking b.c. fllj ] and 
are most relevant for the present content. In this Letter 
a thermodynamic system in the Ising universality class is 
considered in the vicinity of the bulk critical point. The 
system is confined between two parallel plane plates of 



area A separated by distance L. We shall consider only 
those slabs where an external symmetry-breaking bound- 
ary field has been applied to both plates, i.e. a field hi 
(respectively, h%) acting on the plate at z = (respec- 
tively, z = L), and assume that fields hi and hi are of 
the same sign, h\hi > 0, corresponding to the so called 
symmetric b.c. 

Thermodynamic CF is defined as a generalized 
force conjugate to separation L between the plates 
F Casimil (T-L) := where f x (T;L) is the re- 

duced incremental free enegy defined by f x (T;L) := 
lim^oo k B T c A ~ Lfb for free energy F with / b being 
the reduced bulk free energy. It is characterized by the 
property -Fbasimir(r; L) — > as L — + oo. According to 
the finite-size scaling theory, critical phenomena near the 
bulk critical temperature T c and bulk field h = are gov- 
erned by universal scalin g fu nctions that depend on the 
ratio Ll£ [12, 0, ES El, where £ is the bulk corre- 
lation length with £(t,h = 0) rs £q as the reduced 
temperature t = (T — T c )/T c — > , £ nonuniversal 
amplitudes and v a critical exponent. Then the CF can 
be expressed in terms of the universal scaling function 
W ab (-) 1: 

r(T;L) = L~ d * W a b(y)i y = c x tl>' v , 



Casimir I, 



(i) 

where ci is a nonuniversal metric factor. The scaling 
function W a b(y), having universal shape [||, does depend 
on the definition of the correlation length. In order to al- 
low for the "natural" scaling variable L/£ (£ is chosen as 
true correlation length) to emerge in the local-functional 
expressions of W a b{y) m the asymptotic limits y ±oo, 
considered shortly, we choose c\ = l/C^o") 1 ^- Exactly at 
the critical temperature T c , the scaling functions W a b(-) 
give the universal Casimir amplitudes P, Q A ab via 
W a b(0) = (d* — l)A a b, as already considered within local- 
functional theory for symmetric and antisymmetric (H — ) 
(hih 2 < 0) b.c. [4]. Note that d* — min(d, d>), where d is 
a spatial dimension and d> the upper critical dimension 
of the system, the Ising universality class has d > = 4. 
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The purpose of this Letter is t o ap ply EdGF method 
introduced by Fisher and Upton [17|, in order to exam- 
ine the Casimir effect for systems of the Ising universal- 
ity class under the symmetric (H — h) b.c. over the whole 
temperature range, in particular, in d = 3, important 
in respect to the experiments where more accurate theo- 
retical analysis of the above quantities has been missing 
until now. As a nonperturbative approach, EdGF theory 
allows for calculation directly at a fixed spatial dimen- 
sion, an advantage over field-theoretic approach in terms 
of e expansion. 

The local-functional method [I?) asserts that magneti- 
zation profile m(z) in film geometry is given by minimiz- 
ing a (local) interfacial functional JF[m]: 

T[m] := / A{m, m,t, h)dz + f\{mx; hi) + f 2 {m%; h 2 ) 
Jo 

(2) 

where mi = m{z = 0), m, 2 = m(z — L) with f t = 
—hivtii — gmf /2 (i — 1, 2), the usual surface terms which 
allow for the presence of external walls (at z = and 
z = L), and m = dm/dz. The integrand A is assumed to 
take the form which contains only bulk quantities [l7| : 
A(m,m;t,h) = {J(m)Q[A(m,t,h)m] + l}W(m,t,h), 
where W(m,t, h) — <f>(m,t) — $(mb,t) — h(m — f»b), 
and $(m, t) is the bulk Helmholtz free energy density. 
Spontaneous magnetization is denoted by TOb, where 
TOb = B(—t)P for t < and mh = for t > with 
/3 a critical exponent and B a nonuniversal amplitude. 
The function Q(x) is required to satisfy several prop- 
erties 0, Il7| . As before [l7j . we choose J{m) = 1 
and A(m;t,h) := £(m; t)/ \j2x{m; t)W(m; t, h), where 
£(m;i) and x(m;i) are, respectively, the bulk correla- 
tion length and susceptibility of a homogenous system 
at (m,t). Mean-field theory (d > 4) follows from having 
$(m;i) take Landau form with (£ 2 /2x)(m; t) being con- 
stant in m and t. For more general d > 1, bulk functions 
have the following analytic scaling forms [l7j ]: 



W(m;t,0) w |m| a+1 F±(m/TO (t)), 
(C 2 /2x)(m;t) « |m|-^Z±(m/mo(t)), 



(3a) 
(3b) 



in the simultaneous scaling limits i — > and m — ► 0, 
where mo(t) := B\t\P, r\ is the critical bulk correlation 
function exponent in standard notation. 

Minimization of the functional, Eq. ([2]), yields the mag- 
netization profile, m(z), which for hi > 0, contains a 
minimum at z = z+ with magnetization m+ := m(z+) := 
mQ(t)w. The scaling variable y and w are related solely 
in terms of universal quantities: 



Mv\ v 



Z±(u)/Y±(u)du 



u i+"/^|g-i[i - (w/^i+^KtH/Kt^) 



(4) 

where G(x) := xdG/dx-G, A 2 := R x 5/[Q 2 yj28{8 + 1)] is 
defined by the standard universal amplitudes 18] R x — 
C + B s - l D, Q 2 = {C+/C C ){U$?~^ C+ a nonuniver- 
sal zero-field susceptibility amplitude above the critical 



temperature, D, C c and £ c defined along the critical 
isotherm: (T = T c ), h w Dm s , with C c being the corre- 
sponding susceptibility amplitude x ~ C c |/i| - ^ 1-1 /' 5 ' and 
£ c is defined from £(m;0) ~ Cc|^| ^Z' 3 ' 5 . The universal 
functions Y(-) and Z(-) are obtained from normalizing 
Y(-) and Z(-), respectively. The local-functional calcu- 
lation of the CF then follows from df x /dL = W(m+), 
which yields the universal scaling function W++(y) for 
d < 4: 



W ++ (y) = -Ai\y\ 



..1+5 



Y±(w), 



(5) 



with another universal constant A\ := R x (R^) d /[(l + 
S)R C ], defined by other standard universal amplitudes 
[l| R^ = {aA + )^ d ^+ and R c = aA + C+/B 2 , where 
a is the specific-heat exponent and A + a nonuniversal 
specific-heat amplitude for t > and /i = 0. Eqs. (|4|) 
and ([5|) determine completely universal W++(y) within 
the local-functional approach. 

Asymptotic behavior of W ++ (y) as y — > oo follows 
from Eqs. (dJ) and (JSJ) by taking w — > from which Eq. (|U) 
yields y v w 2 ln(B + /iy)+0(i/;), where -B+ is some univer- 
sal constant. Similarly, W++(y) as y — > — oo is obtained 
from (j4|) and ((5|) by taking w — > 1 yielding a similar ex- 
pression for Ij/I" but in terms of w — 1. Solving for w and 
substituting into Eq. ([5]) gives the following 



W++{y) 



-W+^y 2 a exp(-y"), as y 

\2-a 



-oo; 

-^_, 0O |y| 2 - Q exp(-C/ c | 2/ r), asy^-co; 

(6) 

where ?7g = £q /£o , and VK± j0 o are new universal ampli- 
tudes. The results summarized by Eq. are general 
in that that they hold in arbitrary spatial dimension d. 
Previous results, referring to some special cases, such as 
exact calculations on the Ising strip [7| , and on the Ising 
chain subject to two identical surface fields, mean-field 
analysis based on the Ginzburg-Landau tp 4 Hamiltonian 
Q , as well as mean- field treatments of confined fluids [|[ , 
confirm the power-law-exponential behavior of W++(y) 
shown in Eq. ©. 

In obtaining these results one can, to a very good ap- 
proximation, set G(x) — x 2 [4]. This also applies to all 
subsequent results pertaining to the symmetric b.c. and 
greatly simplifies the calculations. 

Mean-field form of W++(y) in terms of the Jacobi 
functions 5] follows also within local-functional approach 
from Eqs. 0$ and (JSJ), when classical values for critical 
exponents are employed along with the scaling functions 
Kfc(-) and Z±{-) for d > 4. 

Excess (Gibbs) adsorption T(t,h). The Gibbs adsorp- 
tion, defined by T(t,h) = J Q [m(z;t, h) — mb(t, h)]dz is 
an integrated measure of the degree of ordering of spins 
or, equivalcntly, in the language of fluids, the amount of 
adsorbed substance on the walls [l!|. From the above 
definition and the scaling postulate H m(z,L,T) w 
mQ(t) r tl> ++ (x, y), x := z/L, valid in the scaling limit 



3 



0, L — > oo, z — > oo, L — z — > do, follows: 
r(i,0) = i?(e + ) /3/ ^ 1 -^G(y), 



G(y) 



:= |tf|" 



(7a) 
(7b) 



with G(y) universal and O(-) the Heaviside function. 
Asymptotically, G{y) ~ |y|' 3 ~ iy as |y| — > oo, so that G(y) 
vanishes for large y and d < 4, j3 < v. Since G(y) is not 
smooth at y = 0, we perfer to express results in terms of 
the universal quantity G(y), defined for d < 4 by 



G(y)=G(y) + \yfe(-y), 

so that J^mdz = B^+f^L^^Giy). 
tional theory predicts that for d < 4: 



G(y) = (l/A 2 )yP-» 



^^du 

Y±(«) U 



(8) 

Local func- 
(9) 



1 v « ' Y±(u) 



For d > A (0 = v = 1/2) scaling forms given by 
Eqs. (|7al7b[) and ([9]) fail and one has to redefine them to 
encompass a logarithmic correction: in the limit L — > oo 
and t -> 0, r(L,T) w B£+ {K x \nL + G{y)) +T {t), with 
ifx as a universal constant, and Fo(i) a non- universal 
additive backround containing both analytic backround 
terms and singular corrections. Mean-field results, that 
follow from Eq. @ when classical values of critical ex- 
ponents are used along with the scaling forms Y+(y) = 
1 + 2/y 2 , Y_(y) = (1 - 1/y 2 ) 2 , Z±(u) = 1 (one reads 
them off from Eqs. ([3a] I3bj) for d > 4) can be expressed 
in terms of a complete elliptic integral of the first kind: 



G(y) 



-V2ln[2K 2 (k)}, y = 4(2k 2 ~l)K 2 {k), 
1/2 < k 2 < 1; 
-V21n[2(l - k 2 )K 2 {k)], 
y = -4(1 + k 2 )K 2 (k) 7 < fc < 1; 
-V^ln^X 2 ^)], y = -4(1 - 2k 2 )K 2 (k), 
< fc 2 < 1/2. 



(10) 

Mean-field universal scaling function G(y) is shown by 
Fig. [H below, together with the d — 3 result. 

To derive quantitative predictions at d = 3 for W ++ (y) 
and G(y) we need to substitute into Eqs. |(SJ) and 
([9]) specific values for bulk critical exponents along with 
suitable choices for Y±(y) and Z±(y). We represent bulk 
scaling functions using parametric models introduced by 
Schofield (20|. These have been developed further [lEEl 
and are believed to give the best available fits to bulk 
data and, by their very construction, to give scaling func- 
tions satisfying required analyticity properties. For our 
purposes, pertaining to the present physical problem sit- 
uated in a one-phase region, the original "linear" para- 
metric model [HI, H(| was found to suffice 22]. At d = 3 
we take j3 — 0.328 and v = 0.632 (all other exponents 
follow from the scaling relations) and a satisfactory fit 



to the bulk amplitude ratios, being properties of bulk 
scaling functions, is provided by taking b 2 — 1.30 and 
<i2 = 0.28 in the notation of [2l| , in the linear model. 
The universal scaling function W++ (y) that follows from 
our calculations in d = 3 is presented in Fig. [T] together 
with an earlier exact curve in spatial dimensions d = 2 
jfl and d > 4 (mean-field) [5(, which not being universal 
23j is shown in reduced universal form VF ++ (y) /W++(0). 
The d = 3 result confirms qualitatively similar structure 




FIG. 1: Plots of the scaling function W++(y) for the Ising 
universality class in (a) d = 3 as obtained by local-functional 
theory (solid line) and compared with MC result [rD] (open 
circles, rectangles and triangles); exact d = 2 (dashed line) 
0; (b) d > 4, presented in a reduced (universal) form derived 
for the first time in [e|. 

of the CF in d = 3 with the ones observed in all the 
other spatial dimensions. This refers to the negative sign 
of the CF for like b.c, smoothness across the whole in- 
terval of scaling variable y £ R (apart from an |y| 2_a 
singularity at y = 0, which can be shown to be quite 
general) as expected based on the fact that the critical 
point of the film (T C (L) , h c (L)) is located o/f the temper- 
ature axis at a non-zero critical bulk field (24|. It also 
follows from this analysis that the minimum of W++ (y) 
in d = 3 is located above the critical point as in other 
dimensions. Fig. Q] shows that there is striking agree- 
ment between the present calculation of (y) and re- 
cent MC simulation results 11], with quoted value of 
W±+(0) = -0.884 implying that the value of the Casimir 
amplitude Ajr9 = —0.442 is closer to this and the earlier 
local-functional result of —0.42(8) [4| as compared to the 
previous MC study [5( . 

Numerical predictions for G(y) in d — 3, based on the 
EdGF Eqs. Q and © within a parametric representa- 
tion, as well as analytic result in the mean-field limit ac- 
cording to Eq. (fTU]) , are given by Fig. [21 showing smooth 
curves, diverging as \y\" for y — > — oo in accord with the 
general definition of G(y). 

There is also much experimental and theoretical inter- 
est in the Casimir effect for the antisymmetric b.c. (-1 — ) 

with recent MC results presented for W-\ (y) [ll|. In 

this case, complications arise in the application of local- 
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FIG. 2: Universal scaling function G(y) of the excess ad- 
sorption r for the Ising universality class calculated within 
local- functional theory (a) in d=3, by Eqs. Q and Q; (b) in 
d > 4 according to the EdGF analytic solution (|10[) . 



functional methods for two main reasons: (i) the approx- 
imation Q{x) — x 2 no longer holds and one needs to use 
the far more complicated form of Q{x) as introduced in 
4]; (ii) one needs to extend the bulk scaling functions 
Y-('), %-{') into the two-phase region, a somewhat ad 
hoc precedure although possible if one uses trigonometric 
parametric models (instead of the linear model) [13, HH 
giving rise to "nonclassical van der Waals loops" . How- 
ever, this more complicated calculation is possible and 
forms the subject of ongoing research. 



More details will follow in a longer report. We kindly 
thank Prof.S. Dietrich and Dr.O. Vasilyev for making 
their MC data [TlJ available for us that enabled com- 
parisons with the present result of EdGF theory within 

Fig.m 
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